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Engineering  and  Design 

INTRODUCTION  TO  PROBABILITY  AND  RELIABILITY  METHODS 
FOR  USE  IN  GEOTECHNICAL  ENGINEERING 


1.  Purpose 

This  engineer  technical  letter  (ETL)  provides  an 
introduction  to  the  use  of  probabilistic  methods  in 
geotechnical  engineering. 


2.  Applicability 

This  ETL  applies  to  HQUSACE  elements  and  USAGE 
commands  having  responsibility  for  the  design  of  civil 
works  projects. 


FOR  THE  DIRECTOR  OF  CIVIL  WORKS: 


A PP  A  -  References 
APP  B  -  Introduction  to  Probability  and 
Reliability  in  Geotechnical 
Engineering 


3.  References 

See  Appendix  A. 


4.  Discussion 

This  is  the  first  in  a  series  of  ETL's  that  will  provide 
guidance  on  the  use  and  application  of  probability  and 
reliability  methods  of  analyses  for  use  in  the 
assessment  of  existing  levees  for  benefit  determination 
and  the  geotechnical  portion  of  major  rehabilitation 
reports. 


Acting  Chief,  Engineering  Division 
Directorate  of  Civil  Works 
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APPENDIX  B:  INTRODUCTION  TO  PROBABILITY  AND 
RELIABILITY  IN  GEOTECHNICAL  ENGINEERING 


B-1.  Introduction 

a.  The  objective  of  this  ETL  is  to  introduce  some 
basic  elements  of  engineering  reliability  analysis 
applicable  to  geotechnical  structures  for  various  modes 
of  performance.  These  reliability  measures  are 
intended  to  be  sufficiently  consistent  and  suitable  for 
application  to  economic  analysis  of  geotechnical 
structures  of  water  resource  projects.  References  are 
provided  which  should  be  consulted  for  detailed 
discussion  of  the  principles  of  reliability  analyses. 

b.  Traditionally,  evaluations  of  geotechnical 
adequacy  are  expressed  by  safety  factors.  A  safety 
factor  can  be  expressed  as  the  ratio  of  capacity  to 
demand.  The  safety  concept,  however,  has  short¬ 
comings  as  a  measure  of  the  relative  reliability  of 
geotechnical  structures  for  different  performance 
modes.  A  primary  deficiency  is  that  parameters 
(material  properties,  strengths,  loads,  etc.)  must  be 
assigned  single,  precise  values  when  the  appropriate 
values  may  in  fact  be  uncertain.  The  use  of  precisely 
defined  single  values  in  an  analysis  is  known  as  the 
deterministic  approach.  The  safety  factor  using  this 
approach  reflects  the  condition  of  the  feature,  the 
engineer's  judgment,  and  the  degree  of  conservatism 
incorporated  into  the  parameter  values. 

c.  Another  approach,  the  probabiiistic  approach, 
extends  the  safety  factor  concept  to  explicitly  incor¬ 
porate  uncertainty  in  the  parameters.  This  uncertainty 
can  be  quantified  through  statistical  analysis  of  existing 
data  or  judgmentally  assigned.  Even  if  judgmentally 
assigned,  the  probabilistic  results  will  be  more 
meaningful  than  a  deterministic  analysis  because  the 
engineer  provides  a  measure  of  the  uncertainty  of  his 
or  her  judgment  in  each  parameter. 


B-2.  Reliability  Analysis  Principles 

a.  The  probabiiity  of  faiiure. 

(1)  Engineering  reliability  analysis  is  concerned 
with  finding  the  reiiabiiity  R  or  the  probabiiity  of 
faiiure  Pr(f)  of  a  feature,  structure,  or  system.  Asa 
system  is  considered  reliable  unless  it  fails,  the 
reliability  and  probability  of  failure  sum  to  unity: 


R  +  Pr[f)  =  1 
R  =  1  -  Pr[f) 

Pr[f)  =  l  R 

(2)  In  the  engineering  reliability  literature,  the 
term  faiiure  is  used  to  refer  to  any  occurrence  of  an 
adverse  event  under  consideration,  including  simple 
events  such  as  maintenance  items.  To  distinguish 
adverse  but  noncatastrophic  events  (which  may  require 
repairs  and  associated  expenditures)  from  events  of 
catastrophic  failure  (as  used  in  the  dam  safety  context), 
the  term  probabiiity  of  unsatisfactory  performance 
Pr(U)  is  sometimes  used.  An  example  would  be  slope 
stability  where  the  safety  factor  is  below  the  required 
minimum  safety  factor  but  above  1.0.  Thus,  for  this 
case,  reliability  is  defined  as: 

R  =  1  -  Pr{U) 


b.  C ontexts  of  reiiabiiity  anaiysis. 

(1)  Engineering  reliability  analysis  can  be  used  in 
several  general  contexts: 

•  The  estimation  of  the  reliability  of  a  new 
structure  or  system  upon  its  construction  and 
first  loading. 

•  The  estimation  of  the  reliability  of  an  existing 
structure  or  system  upon  a  new  loading. 

•  The  estimation  of  the  probability  of  a  part  or 
system  surviving  for  a  given  lifetime. 

N  ote  that  the  third  context  has  an  associated  time 
interval,  where  as  the  first  two  involve  measures  of  the 
overall  adequacy  of  the  system  in  response  to  a  load 
event. 

(2)  Reliability  for  the  first  two  contexts  can  be 
calculated  using  the  capacity-demand  modei  and 
quantified  by  the  reiiabiiity  index  (3.  I  n  the  capacity- 
demand  model,  uncertainty  in  the  performance  of  the 
structure  or  system  is  taken  to  be  a  function  of  the 
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uncertainty  in  the  values  of  various  parameters  used  in 
calculating  some  measure  of  performance,  such  as  the 
factor  of  safety. 

(3)  In  the  third  context,  reliability  over  a  future 
time  interval  is  calculated  using  parameters  developed 
from  actual  data  on  the  lifetimes  or  frequencies  of 
failure  of  similar  parts  or  systems.  These  are  usually 
taken  to  follow  the  exponential  orWeibull  probability 
distributions.  This  methodology  is  well  established  in 
electrical,  mechanical,  and  aerospace  engineering 
where  parts  and  components  routinely  require  periodic 
replacement.  This  approach  produces  a  hazard 
function  which  defines  the  probability  of  failure  in  any 
time  period.  These  functions  are  used  in  economic 
analysis  of  proposed  geotechnical  improvements.  The 
development  of  hazard  functions  is  not  part  of  this 
ETL. 

(4)  For  reliability  evaluation  of  most  geotechnical 
structures,  in  particular  existing  levees,  the  capacity- 
demand  model  will  be  utilized,  as  the  question  of 
interest  is  the  probability  of  failure  related  to  a  load 
event  rather  than  the  probability  of  failure  within  a 
time  interval. 

c.  Reliability  index.  The  reliability  index  (3  is  a 
measure  of  the  reliability  of  an  engineering  system  that 
reflects  both  the  mechanics  of  the  problem  and  the 
uncertainty  in  the  input  variables.  This  index  was 
developed  by  the  structural  engineering  profession  to 
provide  a  measure  of  comparative  reliability  without 
having  to  assume  or  determine  the  shape  of  the 
probability  distribution  necessary  to  calculate  an  exact 
value  of  the  probability  of  failure.  The  reliability 
index  is  defined  in  terms  of  the  expected  value  and 
standard  deviation  of  the  performance  function,  and 
permits  comparison  of  reliability  among  different 
structures  or  modes  of  performance  without  having  to 
calculate  absolute  probability  values.  Calculating  the 
reliability  index  requires: 

•  A  deterministic  model  (e.g.,  a  slope  stability 
analysis  procedure). 

•  A  performance  function  (e.g.,  the  factor  of 
safety  from  UTEXAS2). 

•  The  expected  values  and  standard  deviations 
of  the  parameters  taken  as  random  variables 
(e.g.,  E[(j)j  and  o^). 


•  A  definition  of  the  limit  state  (e.g..  In  (FS)  = 
0). 

•  A  method  to  estimate  the  expected  value  and 
standard  deviation  of  the  limit  state  given  the 
expected  values  and  standard  deviations  of  the 
parameters  (e.g.,  the  Taylor's  series  or  point 
estimate  methods). 

d.  Accuracy  of  reliability  index. 

(1)  For  rehabilitation  studies  of  geotechnical 
structures,  the  reliability  index  is  used  as  a  "relative 
measure  of  reliability  or  confidence  in  the  ability  of  a 
structure  to  perform  its  function  in  a  satisfactory 
manner." 

(2)  The  analysis  methods  used  to  calculate  the 
reliability  index  should  be  sufficiently  accurate  to  rank 
the  relative  reliability  of  various  structures  and 
components.  Flowever,  reliability  index  values  are  not 
absolute  measures  of  probability.  Structures, 
components,  and  performance  modes  with  higher 
indices  are  considered  more  reliable  than  those  with 
lower  indices.  Experience  analyzing  geotechnical 
structures  will  refine  these  techniques. 


B-3.  The  Capacity-Demand  Model 

a.  In  the  capacity-demand  model,  the  probability 
of  failure  or  unsatisfactory  performance  is  defined  as 
the  probability  that  the  demand  on  a  system  or 
component  exceeds  the  capacity  of  the  system  or 
component.  The  capacity  and  demand  can  be 
combined  into  a  single  function  (the  performance 
function),  and  the  event  that  the  capacity  equals  the 
demand  taken  as  the //m/t state.  The  re//at)/7/ty  R  is 
the  probability  that  the  limit  state  will  not  be  achieved 
or  crossed. 

b.  The  concept  of  the  capacity-demand  model  is 
illustrated  for  slope  stability  analysis  in  Figure  B-1. 
Using  the  expected  value  and  standard  deviation  of  the 
random  variables  c  and  (|)  in  conjunction  with  the 
Taylor’s  series  method  or  the  point  estimate  method, 
the  expected  value  and  standard  deviation  of  the 
factor  of  safety  can  be  calculated.  If  it  is  assumed  that 
the  factor  of  safety  is  lognormally  distributed,  then  the 
natural  log  of  the  factor  of  safety  is  normally  distri¬ 
buted.  The  performance  function  is  taken  as  the  log  of 
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Figure  B-1.  The  capacity-demand  model 


the  factor  of  safety,  and  the  limit  state  is  taken  as  the 
condition  In  (F5)  =  0.  The  probability  of  failure  is 
then  the  shaded  area  corresponding  to  the  condition  In 
(F5)  <  0.  If  it  is  assumed  that  the  distribution  on  In 
(F5)  is  normal,  then  the  probability  of  failure  can  be 
obtained  using  standard  statistical  tables. 

c.  Equivalent  performance  functions  and  limit 
states  can  be  defined  using  other  measures,  such  as  the 
exit  gradient  for  seepage. 

d.  The  probability  of  failure  associated  with  the 
reliability  index  is  a  probability  per  structure;  it  has  no 
time-frequency  basis.  Once  a  structure  is  constructed 
or  loaded  as  modeled,  it  either  performs  satisfactorily 
or  not.  N  evertheless,  the  (3  value  calculated  for  an 
existing  structure  provides  a  rational  comparative 
measure. 


B-4.  Steps  in  a  Reliability  Analysis  Using  the 
Capacity-Demand  Model 

As  suggested  by  Figure  B-1  for  slope  stability,  a 
reliability  analysis  includes  the  following  steps: 

•  I  mportant  variables  considered  to  have 
sufficient  inherent  uncertainty  are  taken  as 
random  variables  and  characterized  by  their 
expected  values,  standard  deviations,  and 
correlation  coefficients.  In  concept,  every 


variable  in  an  analysis  can  be  modeled  as  a 
random  variable  as  most  properties  and 
parameters  have  some  inherent  variability  and 
uncertainty.  However,  a  few  specific  random 
variables  will  usually  dominate  the  analysis. 
Including  additional  random  variables  may 
unnecessarily  increase  computational  effort 
without  significantly  improving  results.  When 
in  doubt,  a  few  analyses  with  and  without 
certain  random  variables  will  quickly  illustrate 
which  are  significant,  as  will  the  examination 
of  variance  terms  in  a  T aylor's  series  analysis. 
For  levee  analysis,  significant  random 
variables  typically  include  material  strengths, 
soil  permeability  or  permeability  ratio,  and 
thickness  of  top  stratum.  M  aterial  properties 
such  as  soil  density  may  be  significant,  but 
where  strength  and  density  both  appear  in  an 
analysis,  strength  may  dominate.  An  example 
of  a  variable  that  can  be  represented 
deterministically  (nonrandom)  is  the  density  of 
water. 

•  A  performance  function  and  limit  state  are 
identified. 

•  The  expected  value  and  standard  deviation  of 
the  performance  function  are  calculated.  I  n 
concept,  this  involves  integrating  the  perfor¬ 
mance  function  over  the  probability  density 
functions  of  the  random  variables.  I  n  prac¬ 
tice,  approximate  values  are  obtained  using 
the  expected  value,  standard  deviation,  and 
correlation  coefficients  of  the  random  vari¬ 
ables  in  the  T aylor's  series  method  or  the 
point  estimate  method. 

•  The  reliability  index  (3  is  calculated  from  the 
expected  and  standard  deviation  of  the  perfor¬ 
mance  function.  The  reliability  index  is  a 
measure  of  the  distance  between  the  expected 
value  of  In  (C/D)  or  In  (FS)  and  the  limit 
state. 

•  If  a  probability  of  failure  value  is  desired,  a 
distribution  is  assumed  and  Pr(f)  calculated. 


B-5.  Random  Variables 

a.  Description.  Parameters  having  significance 
in  the  analysis  and  some  significant  uncertainty  are 
taken  as  random  variables.  Instead  of  having  precise 
single  values,  random  variables  assume  a  range  of 
values  in  accordance  with  a  probability  density 
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function  or  probability  distribution.  The  probability 
distribution  quantifies  the  likelihood  that  its  value  lies 
in  any  given  interval.  T wo  commonly  used  distri¬ 
butions,  the  normal  and  the  lognormal,  are  described 
later  in  this  appendix. 

b.  M  oments  of  random  variables.  To  model 
random  variables  in  the  Taylor's  series  or  point 
estimate  methods,  one  must  provide  values  of  their 
expected  values  and  standard  deviations,  which  are 
two  of  several  probabilistic  moments  of  a  random 
variable.  These  can  be  calculated  from  data  or  esti¬ 
mated  from  experience.  For  random  variables  which 
are  not  independent  of  each  other,  but  tend  to  vary 
together,  correlation  coefficients  must  also  be 
assigned. 

(1)  Mean  value.  The  mean  value  of  9  set  of  W 
measured  values  for  the  random  variable  X  is  obtained 
by  summing  the  values  and  dividing  by  W: 


N 

E 

/=i 


(2)  Expected  value.  The  expected  value  E[X]  of  a 
random  variable  is  the  mean  value  one  would  obtain  if 
all  possible  values  of  the  random  variable  were  multi¬ 
plied  by  their  likelihood  of  occurrence  and  summed. 
Where  a  mean  value  can  be  calculated  from 
representative  data,  it  provides  an  unbiased  estimate  of 
the  expected  value  of  a  parameter;  hence,  the  mean 
and  expected  value  are  numerically  the  same.  The  ex¬ 
pected  value  is  defined  as: 


E[X]  =  fix  =  I Xf{X)  dx  =  E  Xp[X) 
where 

f(X)  =  probability  density  function  of  X  (for 
continuous  random  variables) 
p(X,)  =  probability  of  the  value  X,  (for  discrete 
random  variables) 

(3)  Variance.  The  vanance  l/ar  [X]  of  a  random 
variable X  is  the  expected  value  of  the  squared 
difference  between  the  random  variable  and  its  mean 
value.  Where  actual  data  are  available,  the  variance  of 
the  data  can  be  calculated  by  subtracting  each  value 
from  the  mean,  squaring  the  result,  and  determining 
the  average  of  these  values: 


Var[x]  =  £[(X  -  =  [  [X  -  fixfM  ^X 

^  E  [(^,-  -  f^xf] 

N 

The  summation  form  above  involving  theX,  term 
provides  the  variance  of  a  population  containing 
exactly  W  elements.  U  sually,  a  sample  of  size  W  is 
used  to  obtain  an  estimate  of  the  variance  of  the 
associated  random  variable  which  represents  an  entire 
population  of  items  or  continuum  of  material.  To 
obtain  an  unbiased  estimate  of  the  population  working 
from  a  finite  sample,  the  W  is  replaced  by  W  -  1: 

VarlX]  =  ^  ~  ^->'1 

^  ^  N  -  I 


(4)  Standard  deviation.  T o  express  the  scatter  or 
dispersion  of  a  random  variable  about  its  expected 
value  in  the  same  units  as  the  random  variable  itself, 
the  standard  deviation  is  taken  as  the  square  root  of 
the  variance:  thus: 


=  ^Var[Xl 

(5)  Coefficient  of  variation.  To  provide  a 
convenient  dimensionless  expression  of  the  uncertainty 
inherent  in  a  random  variable,  the  standard  deviation  is 
divided  by  the  expected  value  to  obtain  the  coefficient 
of  variation  which  is  usually  expressed  as  a  percent: 


X  100  % 

£[X] 


The  expected  value,  standard  deviation,  and  coefficient 
of  variation  are  interdependent:  knowing  any  two,  the 
third  is  known.  I  n  practice,  a  convenient  way  to 
estimate  moments  for  parameters  where  little  data  are 
available  is  to  assume  that  the  coefficient  of  variation  is 
similar  to  previously  measured  values  from  other  data 
sets  for  the  same  parameter. 


c.  Correlation.  Pairs  of  random  variables  may 
be  correlated  or  independent;  if  correlated,  the  like¬ 
lihood  of  a  certain  value  of  the  random  variable  Y 
depends  on  the  value  of  the  random  variable  X.  For 
example,  the  strength  of  sand  may  be  correlated  with 
density  or  the  top  blanket  permeability  may  be  corre¬ 
lated  with  grain  size  of  the  sand.  The  covariance 
Cov[X,Y]  is  analogous  to  the  variance  but  measures 
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the  combined  effect  of  how  two  variables  vary 
together.  The  definition  of  the  covariance  is: 

Cov[x,y]  =  £[(X  -  -  p,)] 

which  is  equivalent  to: 


Cov  [X,Y]  =  1 1  (X  -  |j^)(F  -  liy)  f{X,Y)  dY  dX 

In  the  above  equation,  f(X,Y)  is  the  joint  probability 
density  function  of  the  random  variables  X  and  T.  To 
calculate  the  covariance  from  data,  the  following 
equation  can  be  used: 

Cov[X.Y\  .  1  Y.  P,  -  -  (I,) 


the  normal  or  lognormal  distribution  which  has  the 
expected  value  and  standard  deviation  assumed.  This 
can  easily  be  done  with  spreadsheet  software. 

(2)  Figure  B-1  illustrated  probability  density 
functions  for  the  random  variables  c  and  cj).  A 
probability  density  function  has  the  property  that  for 
any  X,  the  value  of  f(x)  is  proportional  to  the  like¬ 
lihood  of  X.  The  area  under  a  probability  density 
function  is  unity.  The  probability  that  the  random 
variable  X  lies  between  two  values  Xj  and  Xj  is  the 
integral  of  the  probability  density  function  taken 
between  the  two  values.  H  ence: 


Pr[X^  <  X  <  X^)  =  I  f^{X)  dx 


To  provide  a  nondimensional  measure  of  the  degree  of 
correlation  between  X  and  T,  the  correlation  coefficient 
Pxj  is  obtained  by  dividing  the  covariance  by  the 
product  of  the  standard  deviations: 

Cov 


(3)  The  cumulative  distribution  function  CDF  or 
Fx(X)  measures  the  integral  of  the  probability  density 
function  from  minus  infinity  to  X: 

^xiX)  =  f  fxW  dx 


The  correlation  coefficient  may  assume  values  from 
-1.0  to  -F  1.0.  A  value  of  1.0  or  -1.0  indicates  there  is 
perfect  linear  correlation;  given  a  value  of  X,  the  value 
of  Y  is  known  and  hence  is  not  random.  A  value  of 
zero  indicates  no  linear  correlation  between  variables. 

A  positive  value  indicates  the  variables  increase  and 
decrease  together;  a  negative  value  indicates  that  one 
variable  decreases  as  the  other  increases.  Pairs  of 
independent  random  variables  have  zero  correlation 
coefficients. 


B-6.  Probability  Distributions 

a.  Definition. 

(1)  The  terms  probability  distribution  and  proba¬ 
bility  density  function  pdf  or  the  notation  f^fX)  refer  to 
a  function  that  defines  a  continuous  random  variable. 
The  T aylor's  series  and  point  estimate  methods 
described  herein  to  determine  moments  of  performance 
functions  require  only  the  mean  and  standard  deviation 
of  random  variables  and  their  correlation  coefficients; 
knowledge  of  the  form  of  the  probability  density 
function  is  not  necessary.  However,  in  order  to  ensure 
that  estimates  made  for  these  moments  are  reasonable, 
it  is  recommended  that  the  engineer  plot  the  shape  of 


Thus  for  any  value  X,  F^iX)  is  the  probability  that  the 
random  variable  X  is  less  than  the  given  x. 

b.  Estimating  probabilistic  distributions.  A 
suggested  method  to  assign  or  check  assumed 
moments  for  random  variables  is  to: 

•  Assume  trial  values  for  the  expected  value  and 
standard  deviation  and  take  the  random 
variable  to  be  normal  or  lognormal. 

•  Plot  the  resulting  density  function  and  tabulate 
and  plot  the  resulting  cumulative  distribution 
function  (spreadsheet  software  is  a  convenient 
way  to  do  this). 

•  Assess  the  reasonableness  of  the  shape  of  the 
pdf  and  the  val ues  of  the  CD F . 

•  Repeat  above  steps  with  successively 
improved  estimates  of  the  expected  value  and 
standard  deviation  until  an  appropriate  pdf  and 
CDF  are  obtained. 

c.  Normal  distribution.  The  normal  or  Gaussian 
distribution  is  the  most  well-known  and  widely 
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assumed  probability  density  function.  It  is  defined  in 
terms  of  the  mean  p;(  and  standard  deviation  as: 


MX)  - 


a  \J2ti 


exp 


2o2 


When  fitting  the  normal  distribution,  the  mean  of  the 
distribution  is  taken  as  the  expected  value  of  the 
random  variable.  The  cumulative  distribution  function 
for  the  normal  distribution  is  not  conveniently  ex¬ 
pressed  in  closed  form  but  is  widely  tabulated  and  can 
be  readily  computed  by  numerical  approximation.  It  is 
a  built-in  function  in  most  spreadsheet  programs. 
Although  the  normal  distribution  has  limits  of  plus  and 
minus  infinity,  values  more  than  three  or  four  standard 
deviations  from  the  mean  have  very  low  probability. 
Hence,  one  empirical  fitting  method  is  to  take 
minimum  and  maximum  reasonable  values  to  be  at 
plus  and  minus  three  or  so  standard  deviations.  The 
normal  distribution  is  commonly  assumed  to 
characterize  many  random  variables  where  the 
coefficient  of  variation  is  less  than  about  30  percent. 
For  levees,  these  include  soil  density  and  drained 
friction  angle.  W  here  the  mean  and  standard  deviation 
are  the  only  information  known,  it  can  be  shown  that 
the  normal  distribution  is  the  most  unbiased  choice. 

d.  Lognormal  distribution. 

(1)  When  a  random  variable  /  is  lognormally 
distributed,  its  natural  logarithm,  InX,  is  normally 
distributed.  The  lognormal  distribution  has  several 
properties  which  often  favor  its  selection  to  model 
certain  random  variables  in  engineering  analysis: 

•  As  X  is  positive  for  any  value  of  In  X, 
lognormally  distributed  random  variables 
cannot  assume  values  below  zero. 

•  It  often  provides  a  reasonable  shape  in  cases 
where  the  coefficient  of  variation  is  large 

(>  30  percent)  or  the  random  variable  may 
assume  values  over  one  or  more  orders  of 
magnitude. 

•  The  central  limit  theorem  implies  that  the 
distribution  of  products  or  ratios  of  random 
variables  approaches  the  lognormal  distri¬ 
bution  as  the  number  of  random  variables 
increases. 

(2)  If  the  random  variable  /  is  lognormally 
distributed,  then  the  random  variable  Y  =  In  X  is 
normally  distributed  with  parameters  Ef)']  =  E[ln  X] 


and  oy  =  oinx-  T o  obtain  the  parameters  of  the  normal 
random  variable  Y,  first  the  coefficient  of  variation  of 
X  is  calculated: 


V  =  _ ^ 

^  e[x] 

The  standard  deviation  of  is  then  calculated  as: 

Oy  = 

The  standard  deviation  oy  is  in  turn  used  to  calculate 
the  expected  value  of  Y: 


E[y]  =  £[lnX]  =  ln£[x]  - 


The  density  function  of  the  lognormal  variate  X  is: 


fix) 


exp 


_  1 

'  hi  X  -  E[Yf  ' 

2 

The  shape  of  the  distribution  can  be  plotted  from  the 
above  equation.  Values  on  the  cumulative  distribution 
function  for  X  can  be  determined  from  the  cumulative 
distribution  function  of  Y  (E[Y],  Oy)  by  substituting  the 
X  in  the  expression  Y  =  InX. 


B-7.  Calculation  of  the  Reliability  Index 

a.  Figure  B-2  illustrates  that  a  simple  definition 
of  the  reliability  index  is  based  on  the  assumption  that 
capacity  and  demand  are  normally  distributed  and  the 
limit  state  is  the  event  that  their  difference,  the  safety 
margin  S,  is  zero.  The  random  variable  S  is  then  also 
normally  distributed  and  the  reliability  index  is  the 
distance  by  which  E[SJ  exceeds  zero  in  units  of  0^: 

p  .  £[£]  .  E[c  -D\ 


An  alternative  formulation  (also  shown  in  Figure  B-2) 
implies  that  capacity  C  and  demand  D  are  lognormally 
distributed  random  variables.  In  this  case  In  C  and  In 
D  are  normally  distributed.  Defining  the  factor  of 
safety  FS  as  the  ratio  C/D,  then  In  FS=  (InC)-(ln 
D)  and  In  FS  is  normally  distributed.  Defining  the 
reliability  index  as  the  distance  by  which  In  FS 
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Figure  B-2.  Alternative  definitions  of  the  reliability  index 

exceeds  zero  in  terms  of  the  standard  deviation  of 
In  FS,  it  is: 

£[ln  C  -  InD]  _  £[ln  (C/D)] 


^(In  C  -  In  D  ) 

_  £[lnF5] 


®  In  (C/D) 


'Inf 5 


b.  F rom  the  properties  of  the  lognormal 
distribution,  the  expected  value  of  In  C  is: 


£[lnC]  =  ln£[C]  -  lal. 


where: 


OlnC  =  In 


1  +  v;: 


Similar  expressions  apply  to  E[ln  D]  and  a|„o. 

c.  The  expected  value  of  the  log  of  the  factor  of 
safety  is  then: 


£[lnF5]  =  ln£[C]  -  ]n  E[d]  “  ^  In 


1  +  V" 


+  -  In 
2 


1  +  V" 


As  the  second-order  terms  are  small  when  the  coeffi¬ 
cients  of  variation  are  not  exceedingly  large  (below 
approximately  30  percent),  the  equation  above  is 
sometimes  approximated  as: 


E[\nFS]  =  \nE[C]  -  In  E[d]  =  In  (  |||| 

The  standard  deviation  of  the  log  of  the  factor  of 
safety  is  obtained  as: 


'lnf5 


f 


InC  ^  ^\nD 


'Inf 5 


f 


=  \/ln 


1  +  v; 


+  In 


1  + 


Introducing  an  approximation. 


In 


1  +  Vr 


the  reliability  index  for  lognormally  distributed  C,  D, 
and  FS  and  normally  distributed  In  C,  InD,  and  In  FS 
can  be  expressed  approximately  as: 


In 


P 


E[C\ 

E[D\ 


/ 


The  exact  expression  is: 

E\C]f 


P  =  In 


+ 


:[D]f 


+  V, 


C  J 


For  many  geotechnical  problems  and  related  deter¬ 
ministic  computer  programs,  the  output  is  in  the  form 
of  the  factor  of  safety,  and  the  capacity  and  demand 
are  not  explicitly  separated.  The  reliability  index  must 
be  calculated  from  values  of  E[FS]  and  Ofs  obtained 
from  multiple  runs  as  described  in  the  next  section.  In 
this  case,  the  reliability  index  is  obtained  using  the 
following  steps: 

o. 


V  = 


"FS 


£[lnF5] 


'lnf5 


1  +  V, 


FS 


ln£[F5]  -  |ln(l  +  4) 
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p  _  Eflii  FS] 

^  In  FS 


In  £ 


:[re]  /  / 


1  *  V, 


FS 


In  1  +  V, 


FS 


B-8.  Integration  of  the  Performance  Function 

M  ethods  such  as  direct  integration,  T ayior's  series 
method,  point  estimate  method,  and  M  onte  Cario 
simuiation  are  avaiiabie  for  caicuiating  the  mean  and 
standard  deviation  of  the  performance  function.  F or 
direct  integration,  the  mean  vaiue  of  the  function  is 
obtained  by  integrating  over  the  probabiiity  density 
function  of  the  random  variabies.  A  brief  description 
of  the  other  methods  foiiows.  The  references  shouid 
be  consuited  for  additionai  information. 

a.  Taylor's  series  method.  Tayior's  series 
method  is  one  of  severai  methods  to  estimate  the 
moments  of  a  performance  function  based  on  moments 
of  the  input  random  variabies.  It  is  based  on  a  Taylor's 
series  expansion  of  the  performance  function  about 
some  point.  FortheCorps'  navigation  rehabilitation 
studies,  the  expansion  is  performed  about  the  expected 
values  of  the  random  variables.  The  Taylor's  series 
method  is  termed  a  first-order,  second-moment 
(FOSM  )  method  as  only  first-order  (linear)  terms  of 
the  series  are  retained  and  only  the  first  two  moments 
(mean  and  the  standard  deviation)  are  considered.  The 
method  is  summarized  below. 

(1)  Independent  random  variables.  Given  a 
function  Y  =  g(Xp  Xj,  ...  X„j,  where  all  X,  values  are 
independent,  the  expected  value  of  the  function  is 
obtained  by  evaluating  the  function  at  the  expected 
values  of  the  random  variables: 

E[Y]  =g(£[Xi],  £[x4  ...£[X„]) 

For  a  function  such  as  the  factor  of  safety,  this  implies  that 
the  expected  value  of  the  factor  of  safety  is  calculated 
using  the  expected  values  of  the  random  variables: 


=  F5(£[(i)i],  £[ci],  £[Yi],  - ) 

The  variance  of  the  performance  function  is  taken  as: 


Var[Y]  =  Y. 


dY 

dx, 


*  / 


VarX. 


with  the  partial  derivatives  taken  at  the  expansion  point 
(in  this  case  the  mean  or  expected  value).  U  sing  the 
factor  of  safety  as  an  example  performance  function, 
the  variance  is  obtained  by  finding  the  partial  deriva¬ 
tive  of  the  factor  of  safety  with  respect  to  each  random 
variable  evaluated  at  the  expected  value  of  that  varia¬ 
ble,  squaring  it,  multiplying  it  by  the  variance  of  that 
random  variable,  and  summing  these  terms  over  all  of 
the  random  variables: 


Var[FS]  =  Y 


The  standard  deviation  of  the  factor  of  safety  is  then 
simply  the  square  root  of  the  variance. 

(a)  FI  aving  the  expected  value  and  variance  of  the 
factor  of  safety,  the  reliability  index  can  be  calculated 
as  described  earlier  in  this  appendix.  Advantages  of 
the  T ayior's  series  method  include  the  following: 

•  The  relative  magnitudes  of  the  terms  in  the 
above  summation  provide  an  explicit  indica¬ 
tion  of  the  relative  contribution  of  uncertainty 
of  each  variable. 

•  The  method  is  exact  for  linear  performance 
functions. 

Disadvantages  of  the  T ayior's  series  method  include 
the  following: 

•  It  is  necessary  to  determine  the  value  of 
derivatives. 

•  The  neglect  of  higher  order  terms  introduces 
errors  for  nonlinear  functions. 

(b)  The  required  derivatives  can  be  estimated 
numerically  by  evaluating  the  performance  function  at 
two  points.  The  function  is  evaluated  at  one  increment 
above  and  below  the  expected  value  of  the  random 
variable  X,  and  the  difference  of  the  results  is  divided 
by  the  difference  between  the  two  values  of  X,. 
Although  the  derivative  at  a  point  is  most  precisely 
evaluated  using  a  very  small  increment,  evaluating  the 
derivative  over  a  range  of  plus  and  minus  one  standard 
deviation  may  better  capture  some  of  the  nonlinear 
behavior  of  the  function  over  a  range  of  likely  values. 
Thus,  the  derivative  is  evaluated  using  the  following 
approximation: 


dFS 

dx, 


VarX. 
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dx.  2o„ 

I  X. 

When  the  above  expression  is  squared  and  multiplied 
by  the  variance,  the  standard  deviation  term  in  the 
denominator  cancels  the  variance,  leading  to 


that  of  the  continuous  random  variable.  Having  only  a 
few  values  over  which  to  integrate,  the  moments  of  the 
performance  function  are  easily  obtained.  A  simple 
and  straightforward  point  estimate  method  has  been 
proposed  by  Rosenblueth  (1975,  1981)  and  is 
summarized  by  Harr  (1987).  That  method  is  briefly 
summarized  below. 


/ 


V 


by" 

dX. 

^  ) 


VarX  = 


kp.)  -  «p->r 
2 


where  X^.  and  X.  are  values  of  the  random  variable  at 
plus  and  minus  one  standard  deviation  from  the 
expected  value. 

(2)  Correlated  random  variables.  Where  random 
variables  are  correlated,  the  solution  is  more  complex. 
The  expression  for  the  expected  value,  retaining 
second-order  terms  is: 


(1)  Independent  random  variables.  As  shown  in 
Figure  B-3,  a  continuous  random  variable X  is 
represented  by  two  point  estimates,  X+  and  X.,  with 
probability  concentrations  P+  and  P.,  respectively.  As 
the  two  point  estimates  and  their  probability  concen¬ 
trations  form  an  equivalent  probability  distribution  for 
the  random  variable,  the  two  P  values  must  sum  to 
unity.  The  two  point  estimates  and  probability  con¬ 
centrations  are  chosen  to  match  three  moments  of  the 
random  variable.  When  these  conditions  are  satisfied 
for  symmetrically  distributed  random  variables,  the 
point  estimates  are  taken  at  the  mean  plus  or  minus  one 
standard  deviation: 


E[Y]  =  g 


:[x,l  e\xA...e\x  ] 

IJ’  [  2J’  [  n 


+  -  E  ^  ^  CovlX.X  ) 
2  dXAX.  ^ 

’  J 


However,  in  keeping  with  the  first-order  approach,  the 
second-order  terms  are  generally  neglected,  and  the 
expected  value  is  calculated  the  same  as  for 
independent  random  variables.  The  variance, 
however,  is  taken  as: 


Var[Y] 


E 


dY 

dx, 


2 

VarX. 


-2E 


dY  dY 
dX.  dX. 

'  J 


c».(x,x,) 


where  the  covariance  term  contains  terms  for  each 
possible  combination  of  random  variables. 

b.  Point  estimate  method.  An  alternative  method 
to  estimate  moments  of  a  performance  function  based 
on  moments  of  the  random  variables  is  the  point 
estimate  method.  Point  estimate  methods  are  pro¬ 
cedures  where  probability  distributions  for  continuous 
random  variables  are  modeled  by  discrete  "equivalent" 
distributions  having  two  or  more  values.  The  elements 
of  these  discrete  distributions  (or  point  estimates)  have 
specific  values  with  defined  probabilities  such  that  the 
first  few  moments  of  the  discrete  distribution  match 


X,-  ■  Ep,]  - 


Figure  B-3.  Point  estimate  method 

For  independent  random  variables,  the  associated 
probability  concentrations  are  each  one-half: 

Pi.  =  Pi-  =  0-50 

Knowing  the  point  estimates  and  their  probability 
concentrations  for  each  variable,  the  expected  value  of 
a  function  of  the  random  variables  raised  to  any 
power  M  can  be  approximated  by  evaluating  the 
function  for  each  possible  combination  of  the  point 
estimates  (e.g.,  X^, ,  Xj. ,  X3+ ,  X„. ),  multiplying  each 
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result  by  the  product  of  the  associated  probability 
concentrations  (e.g.,  P  +  ..  =  Pi+  P2-P3. )  and  summing 
the  terms.  For  example,  two  random  variables  result 
in  four  combinations  of  point  estimates  and  four  terms: 

e[y^]  =  +  P^_g[X,^,X^_Y 

For  W  random  variables,  there  are  2*“  combinations  of 
the  point  estimates  and  2*“  terms  in  the  summation.  To 
obtain  the  expected  value  of  the  performance  function, 
the  function  g(XyX2)  is  calculated  2*“  times  using  all  the 
combinations  and  the  exponent  M  is  1.  To  obtain  the 
standard  deviation  of  the  performance  function,  the 
exponent  M  is  taken  as  2  and  the  squares  of  the 
obtained  results  are  weighted  and  summed  to  obtain 
E[Y^].  The  variance  can  then  be  obtained  from  the 
identity 

Var[Y]  =  e[y^]  - 

and  the  standard  deviation  is  the  square  root  of  the 
variance. 

(2)  Correlated  random  variables.  Correlation 
between  symmetrically  distributed  random  variables  is 
treated  by  adjusting  the  probability  concentrations 
(P  ±  ±  ....  ±  ).  A  detailed  discussion  is  provided  by 
Rosenblueth  (1975)  and  summarized  by  Flarr  (1987). 
For  certain  geotechnical  analyses  involving  lateral 
earth  pressure,  bearing  capacity  of  shallow  founda¬ 
tions,  and  slope  stability,  often  only  two  random 
variables  (c  and  cj)  or  tan  cj))  need  to  considered  as 
correlated.  For  two  correlated  random  variables 
within  a  group  of  two  or  more,  the  product  of  their 
concentrations  is  modified  by  adding  a  correlation 
term: 


B-9.  Determining  the  Probability  of  Failure 

Once  the  expected  value  and  standard  deviation  of  the 
performance  function  have  been  determined  using  the 
Taylor’s  series  or  point  estimate  methods,  the  reli¬ 
ability  index  can  be  calculated  as  previously  described. 
If  the  reliability  index  is  assumed  to  be  the  number  of 
standard  deviations  by  which  the  expected  value  of  a 
normally  distributed  performance  function  (e.g.. 

In  (PS))  exceeds  zero,  than  the  probability  of  failure 
can  be  calculated  as: 

Pr{f)  =  il;(-p)  =  il;(-z) 

where  ijjf-z)  is  the  cumulative  distribution  function  of 
the  standard  normal  distribution  evaluated  at-z,  which 
is  widely  tabulated  and  available  as  a  built-in  function 
on  modern  microcomputer  spreadsheet  programs. 


B-10.  Overall  System  Reliability 

Reliability  indices  for  a  number  of  components  or  a 
number  of  modes  of  performance  may  be  used  to 
estimate  the  overall  reliability  of  an  embankment. 
There  are  two  types  of  systems  that  bound  the  possible 
cases,  the  series  system  and  the  parallel  system. 

a.  Series  system.  In  a  series  system,  the  system 
will  perfrom  unsatisfactorily  if  any  one  component 
performs  unsatisfactorily.  If  a  system  has  n  compo¬ 
nents  in  series,  the  probability  of  unsatisfactory 
performance  of  the  ith  component  is  p,  and  its  reli¬ 
ability,  R,  =  1  -  p,,  then  the  reliability  of  the  system, 
or  probability  that  all  components  will  perform  satis¬ 
factorily,  is  the  product  of  the  component  reliabilities. 


p.  . 

=  p.  = 

(p.  \(p.  \ 

_  p 

l+J- 

1-J+ 

V  *-/v  ;+/ 

4 

p 

=  p.  = 

IP.  \ 

+  P 

*+J+ 

l-J- 

V  '  +  /v  J+/ 

4 

c.  M onte  Carlo  simulation.  The  performance 
function  is  evaluated  for  many  possible  values  of  the 
random  variables.  A  plot  of  the  results  will  produce 
an  approximation  of  the  probability  distribution.  Once 
the  probability  distribution  is  determined  in  this 
manner,  the  mean  and  standard  deviation  of  the 
distribution  can  be  calculated. 


=  (1  -  /'iHi  -  Pip  -  -  P>}-i^  ^  P: 


b.  Simple  parallel  system.  In  a  parallel  system, 
the  system  will  only  perform  unsatisfactorily  if  all 
components  perform  unsatisfactorily.  Thus,  the 
reliability  is  unity  minus  the  probability  that  all  com¬ 
ponents  perform  unsatisfactorily,  or 

P  =  1  p^p^p^...p....p^ 
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c.  Parallel  series  systems. 


B-11.  Target  Reliability  Indices 


(1)  Solutions  are  available  for  systems  requiring 
r-out-of-n  operable  components,  which  may  be 
applicable  to  problems  such  as  dewatering  with 
multiple  pumps,  where  r  is  defined  as  the  number  of 
reliable  units.  Subsystems  involving  independent 
parallel  and  series  systems  can  be  mathematically 
combined  by  standard  techniques. 

(2)  U  pper  and  lower  bounds  on  system  reliability 
can  be  determined  by  considering  all  components  to  be 
from  subgroups  of  parallel  and  series  systems,  respec¬ 
tively:  however,  the  resulting  bounds  may  be  so  broad 
as  to  be  impractical.  A  number  of  procedures  are 
found  in  the  references  to  narrow  the  bounds. 

(3)  Engineering  systems  such  as  embankments  are 
complex  and  have  many  performance  modes.  Some  of 
these  modes  may  not  be  independent;  for  instance, 
several  performance  modes  may  be  correlated  to  the 
occurrence  of  a  high  or  low  pool  level.  Rational 
estimation  of  the  overall  reliability  of  an  embankment 
is  a  topic  that  is  beyond  this  ETL. 

d.  A  practical  approach. 

(1)  The  reliability  of  a  few  subsystems  or  com¬ 
ponents  may  govern  the  reliability  of  the  entire  system. 
Thus,  developing  a  means  to  characterize  and  compare 
the  reliability  of  these  components  as  a  function  of 
time  is  sufficient  to  make  engineering  judgments  to 
prioritize  operations  and  maintenance  expenditures. 

(2)  For  initial  use  in  reliability  assessment  of 
geotechnical  systems,  the  target  reliability  values 
presented  below  should  be  used.  The  objective  of  the 
M  ajor  Rehabilitation  Program  would  be  to  keep  the 
reliability  index  for  each  significant  mode  above  the 
target  value  for  the  foreseeable  future. 


Reliability  indices  are  a  relative  measure  of  the  current 
condition  and  provide  a  qualitative  estimate  of  the 
expected  performance.  Embankments  with  relatively 
high  reliability  indices  will  be  expected  to  perform 
their  function  well.  Embankments  with  low  reliability 
indices  will  be  expected  to  perform  poorly  and  present 
major  rehabilitation  problems.  If  the  reliability  indices 
are  very  low,  the  embankment  may  be  classified  as  a 
hazard.  The  target  reliability  values  shown  in 
Table  B-1  should  be  used  in  general. 


Table  B-1. 

Target  Reliability  Indices 

Expected 

Performance 

Level 

Beta 

Probability  of 

Unsatisfactory 

Performance 

High 

5 

0.0000003 

Good 

4 

0.00003 

Above  average 

3 

0.001 

Below  average 

2.5 

0.006 

Poor 

2.0 

0.023 

Unsatisfactory 

1.5 

0.07 

Hazardous 

1.0 

0.16 

Note:  Probability  of  unsatisfactory  performance  is  the  probability 
that  the  value  of  performance  function  will  approach  the  limit 
state,  or  that  an  unsatisfactory  event  will  occur.  For  example,  if 
the  performance  function  is  defined  in  terms  of  slope  instability, 
and  the  probability  of  unsatisfactory  performance  is  0.023,  then 
23  of  every  1 ,000  instabilities  will  result  in  damage  which  causes 
a  safety  hazard. _ 
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